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Abstract. Model atmospheres of A and B stars are computed taking into account magnetic line blanketing. These calculations 
are based on the new stellar model atmosphere code LLModels which implements direct treatment of the opacities due to the 
bound-bound transitions and ensures an accurate and detailed description of the line absorption. The anomalous Zeeman effect 
was calculated for the field strengths between 1 and 40 kG and a field vector perpendicular to the line of sight. The model 
structure, high-resolution energy distribution, photometric colors, metallic line spectra and the hydrogen Balmer line profiles 
are computed for magnetic stars with different metallicities and are discussed with respect to those of non-magnetic reference 
models. The magnetically enhanced line blanketing changes the atmospheric structure and leads to a redistribution of energy 
in the stellar spectrum. The most noticeable feature in the optical region is the appearance of the 5200 A depression. However, 
this effect is prominent only in cool A stars and disappears for higher effective temperatures. The presence of a magnetic field 
produces opposite variation of the flux distribution in the optical and UV region. A deficiency of the UV flux is found for the 
whole range of considered effective temperatures, whereas the "null wavelength" where flux remains unchanged shifts towards 
the shorter wavelengths for higher temperatures. 
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1. Introduction 

Magnetic chemically peculiar (CP) stars are the upper and mid- 
dle main sequence stars characterized by anomalous chemi- 
cal abundances and an unusual distribution of energy in their 
spectra. r ei f of magnetic CP stars ranges from 6500 K to about 
25 000 K, whereas the typical magnetic field strengths observed 
in their atmospheres span from a few hundred gauss to ~ 30 kG. 
A strong magnetic field is expected to have an important effect 
on the atmospheres and spectra of CP stars. Whereas s everal 
detailed line profile studies (e.g. Kochukh ov et alJ l2004) have 
included the full treatment of the Zeeman effect in calculation 
of the Stokes parameters of individual spectral lines and short 
wavelength regions, an overall influence of the magnetic ab- 
sorption on the structure of the model atmospheres and flux 
distributions of CP stars remained relatively unexplored. 

Several characteristic features in the energy distribution of 
magnetic CP stars are suspected to be a result of the enhanced 
line blanketing due to the m agnetic intensification of spectral 
lines. For instance, Kodaira ( 1969) found several flux degres- 
sions in the visual spectrum of CP stars, and lLeckrond dl973h 
showed that the UV flux of CP stars is depressed compared to 
that of normal stars with a similar flux distribution in the opti- 
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cal region. These possible manifestations of the magnetic line 
blanketing emphasize a necessity to consider a magnetic field 
in the model atmosphere calculation of CP stars. 

In general, a magnetic field influences the energy transport, 
hydrostatic equilibrium, diffusion proc esses and the formation 
of sp e ctral lines. Previ o us studies dStepien| ll978t Muthsam 
| 1979| ICarnenteJ fl985l iLeBlanc et alJ 1 19941 IValvavin et alJ 
2004) have made attempts to model some of these factors. 
However, due to a limitation of computer resources, it was im- 
possible to fully account for the Zeeman effect on the line ab- 
sorption, except in the modelling of the pure hydrogen atmo- 
spheres of magnetic white dwarfs ( Wickramasinghe & Martin 
1986J). In the early model atmosphere calculations for non- 
degenerate stars magnetic splitting was treated very approxi- 
mately by introducing a pseudo-microturb ulent velocity which 
affects all spectral lines in the same way (Muthsam 1979) or by 
adopting a normal Zeeman triplet pattern with the same effec- 
tive Lande factors (Carpenter 1985) for all lines. Recent inves- 
tigation by|Stift & Leone (2003) demonstrated that the mag- 
netic intensification of spectral lines depends primarily on the 
parameters of the anomalous Zeeman splitting pattern, in par- 
ticular on the number of Zeeman components. In the light of 
these results it becomes clear that previous attempts to simu- 
late magnetic line blanketing by an enhanced microturbulence 
or by using a simple triplet pattern are insufficient. 

The primary aim of our paper is to introduce a realistic 
calculation of the anomalous Zeeman effect in the classical 
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1-D models of stellar atmospheres and to investigate the re- 
sulting effects on the model structure, energy distribution and 
other common observables. We calculate a grid of model at- 
mospheres of A and B stars assuming a horizontal magnetic 
field and exploring a T e ff range relevant for CP stars. The pa- 
per is organized as follows. In Sect.[5]we describe our model 
atmosphere code and numerical implementation of the Zeeman 
effect in the line opacity calculation. Sect.|3]presents numerical 
results and Sect.|4]summarizes our work and compares the ob- 
served and predicted anomalies in the shape and variability of 
the flux distribution of magnetic CP stars. 

2. Calculation of magnetic model atmospheres 

2.1. The stellar model atmosphere code LLModels 

In calculations described in the present paper we employed 
the stellar model a tmosphere code LLModels developed by 
IShulvak et al. (2004). This code uses a direct method, the so- 
called line-by-line or LL technique, for the line opacity calcu- 
lation. Such an approach allowed us to account for the anoma- 
lous Zeeman splitting of spectral lines in the line blanketing 
calculation and, hence, to achieve a qualitatively new level of 
accuracy in modelling the influence of magnetic field on the 
stellar atmospheric structure. 

LLModels is the stellar model atmosphere code which 
aims at modelling the early and intermediate type stars tak- 
ing into account their individual chemical composition and 
an inhomogeneous vertical distribution of elemental abun- 
dances. The c ode is based on the modified atlas9 subroutines 
<Kuruczll993h and on the spectrum synthesis code described by 
iTsvm bal ( 1996). The LLModels code is written in Fortran 90 
and uses the following general approximations: 

- the plane-parallel geometry is assumed; 

- the Local Thermodynamic Equilibrium (LTE) is used to 
calculate the atomic level populations for all chemical 
species; 

- the stellar atmosphere is assumed to be in a steady state; 

- the radiative equilibrium condition is fulfilled. 

The main goal of the LL modelling technique is to avoid a 
simplified or statistical description of the bound-bound opac- 
ity in stellar atmospheres. The temperature distribution for a 
given model strongly depends upon the accuracy of the flux in- 
tegration over the whole spectral range where the star radiates 
significantly. Consequently, to ensure an accurate flux integra- 
tion, the number of frequency points must be sufficient to pro- 
vide a realistic description of individual spectral lines. In the 
past limited computer resources and imperfect numerical tech- 
niques were required to solve this problem. The best known 
appro ximations are the Opacity Distribution Function method 
(ODF, Kurucz 1979) and the Opacity Sampling technique (OS, 
iGustafsson et all 1975b . 

The key idea of the ODF method is to replace the line opac- 
ity at a given wavelength interval AA by a smooth distribution 
function f(AA). The ODF tables are calculated for a set of T -P 
pairs, chemical compositions (solar or scaled solar) and micro- 
turbulent velocities. A sparse wavelength grid is sufficient for 



the flux integration and, thus, the method achieves a substantial 
saving of computing time. On the other hand, the OS method 
relies on statistically random distribution of wavelength points 
over the whole spectral range. A sufficient accuracy is consid- 
ered to be obtained when the flux integral becomes independent 
from the choice of the wavelength grid within a given error. 

Due to their statistical nature, both the ODF and the OS 
methods suffer from serious shortcomings. For example, in the 
case of the ODF technique it is necessary to recalculate the 
ODF tables each time when new stellar abundances are used 
or an input line list is modified due to the magnetic effects. 
Obviously, this method is not efficient in the context of the 
present study because of the necessity to compute the ODF ta- 
bles for a large number of magnetic field strengths and chem- 
ical compositions. On the other hand, the accuracy of the OS 
method directly depends on the number of frequency points 
and their distribution over the considered spectral range. This 
requires a careful calibration of the method each time an im- 
portant input parameter, such as the magnetic field strength, is 
modified. 

With the LL method we do not use any precalculated opac- 
ity tables and, hence, do not require an excessive hard disk 
storage space. Furthermore, no approximations about the line 
opacity coefficient are being made during the model calcula- 
tion. By using a large number of wavelength points 300000 
- 500 000) we provide a detailed description of the line absorp- 
tion for all atmospheric depths and achieve a higher dynamical 
range in opacity. This allows to reach a required accuracy of the 
resulting model atmosphere structure, especially in the upper 
atmospheric layers. Due to the efficient numerical algorithms 
implemented in LLModels, the code is able to compute model 
atmospheres with m odern PCs in a reasonable amount of time 
dShulvak et alJ2004 . 

In all calculations presented here we employed two criteria 
to control the convergence of models: the constancy of the total 
flux and conservation of the radiative equilibrium. Both criteria 
are checked after each iteration and for each atmospheric depth. 
If the condition 

(#rad + ffconv) - 0-T* s = (1) 

and the total energy balance 

j (a v + e Y )J v dv - j(a v + ( V )S v dv = (2) 

are satisfied in all layers with errors of less than 1%, we as- 
sume that a model has converged. Here the and H com rep- 
resent the radiative and convective fluxes, T e ff is the effective 
temperature, a v and ( v are the continuum and line absorption 
coefficients at the frequency v, whereas the J v and S v are the 
mean intensity and source functions respectively. We note that 
the first criterion is not very sensitive to the temperature varia- 
tions in the optically thin layers. Is is the radiative equilibrium 
condition which becomes more relevant in this case. If any of 
these criteria cannot be satisfied, either the initial model param- 
eters were inappropriate or the model does not converge due to 
fundamental limitations of the adopted physical description of 
the stellar atmospheric processes. 
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Table 1. Description of the spectral line lists used in the line blanketing calculations for different effective temperature T e ff and 
metallicity [M/H] = lg(N meta i s /NH) - lg{N mela i s /Nu)e - N pies is the number of spectral lines retained by the preselection procedure, 
A^Lande pres - the number of spectral lines with known Lande factors, A^LandeLS - the number of spectral lines of light elements 
for which Lande factors were calculated assuming LS coupling, A^ande triplet - the number of spectral lines of light and heavy 
elements without a proper term designation in VALD which were splitted assuming a classical Zeeman triplet pattern with the 
effective Lande factor 1.2, N tota i - the total number of the Zeeman components included in the line blanketing calculation. 





[M/H] 


N 

1 ' pres 


*V Lande pres 


A^Lande LS 


Lande triplet 


Motel 


8000 


0.0 


333 940 


305 109 


16 968 


11863 


5 910 394 




0.5 


491 340 


455 591 


20371 


15 378 


8 774132 




1.0 


719988 


676 889 


23 695 


19 404 


12983 819 


11000 


0.0 


284 361 


263 739 


14462 


6160 


5 187013 




0.5 


418 093 


390641 


18 973 


8479 


7 707 542 




1.0 


624273 


588 005 


24 285 


11983 


11625 934 


15 000 


0.0 


346 912 


326 327 


15428 


5157 


6438 224 




0.5 


519 902 


493 198 


19114 


7590 


9731721 




1.0 


785 203 


751050 


23 542 


10611 


14 827 741 



The LLModels code uses either atomic line lists compiled 
bvlKuruczH 19931) or the Vienna Atomic Line Database (VALD, 
iKuokaet al.l fl999) line list which is converted to a special bi- 
nary format. In model atmosphere calculations there is no need 
to use all available spectral line data and include very weak 
lines. Hence, the LLModels code relies on a preselection pro- 
cedure to choose spectral lines contributing significantly to the 
total absorption coefficient. The code selects spectral lines for 
which { v /a v > s, where s is the adopted selection threshold 
and l v and a v are defined above. 



2.2. Line lists 

The magnetic line blanketing was taken into account for all 
spectral lines except the hydrogen lines according to the in- 
dividual anomalous Zeeman splitting pattern. The initial line 
lists were extracted from VALD. The total number of spec- 
tral lines, including lines originating from the predicted lev- 
els, was more than 21.6 million for the spectral range between 
50 and 100 000 A. This line list was used for the preselection 
procedure in the LLModels code using the selection thresh- 
old s = 1%. Preselection allowed us to decrease the number 
of spectral lines involved in the line blanketing calculation to 
about 300000-800000 depending on the model atmosphere 
parameters. 

To calculate a Zeeman splitting of a spectral line one has 
to know the Lande factors of the lower and upper atomic lev- 
els of the respective transition. The VALD database does not 
provide Lande factors for all lines. We found that roughly 4- 
10% of the preselected spectral lines lack information about 
the Lande factors. For the light elements, from He to Sc, the 
LS coupling approximation is usually sufficient and was em- 
ployed where necessary using the term designation provided 
for each line in the VALD line lists. This approach reduced 
the number of lines with missing Lande factors to 1-4%. For 
these remaining lines without a proper term designation in the 
VALD lists we assumed a classical Zeeman triplet splitting pat- 
tern with the effective Lande factor g e ff = 1.2 (which is the 
average value for the spectral lines with computed or experi- 



mental Lande factors). Further information on the statistics of 
the considered spectral line lists is presented in Tabled 

2.3. Zeeman effect in the line opacity 

In the presence of a magnetic field an atomic level k defined by 
the quantum numbers J k , L k , S k splits into the 2J k + 1 states 
with the magnetic quantum numbers M k = -J k , . . . , +J k . The 
absolute value of the splitting is defined by the field modulus 
\B\ and by the Lande factor g k , which in the case of LS coupling 
can be calculated as 



8k =2 + 



3 S k (S k + l)-L k (L k +l) 



2J k (J k + 1) 



(3) 



According to the selection rules the following transitions are 
allowed between the splitted upper u and lower I levels 



C +i = b 

AM = M u - M, = I = p. 

-1 =r 



(4) 



For a normal Zeeman triplet the subscript p corresponds to the 
unshifted n component, whereas b and r correspond to the blue- 
and red-shifted <x components respectively. In the general case 
of an anomalous splitting pattern the indices p, b, r refer to the 
groups of the n and cr components. 

The wavelength shift of the component ij (j = p, b, r, 
ip,b,r = 1, • ■ ■ , N p j,y, N p fi,r is the number of components in each 
group) relative to the laboratory line centre Aq is defined by 



Mi. = 



Amnc 1 



(giM, - g u M u ) h . 



(5) 



Th e relative s trengths S ,■ of the n and cr components are given 
bv lSobelmar] ( ll977l) and' are listed in Tabled 

In general, to calculate a stellar spectrum in the presence 
of a magnetic field one has to consider the polarized radiative 
transfer equation for the Stokes IQUV parameters. A solution 
of this problem requires special num erical techniques and is 
very computationally expensive (e.g. Pisk unov & Kochuk hov 
2002) and, hence, cannot be easily included in the routine 
model atmosphere calculation of magnetic stars. Nevertheless, 
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Table 2. Relative strengths of the Zeeman n and <x components. 
Here AJ = J u - Ji and M-, = -J u , . . . , +J U . In the Zeeman 
regime considered in our paper symmetric components have 
the same strength, e.g. S 1 ,; = S i b . 



AJ 


St, 


Si b , 







M' 


i/4(y„ 


- M,)(J U + 1 + Mj) 


+ 1 


4-Mj 


l/4(/„ 


- Mi)(J u - 1 - Mj) 


-1 


(J u + l) 2 - Mf 


i/4(y„ 


+ 1 + M,)(J U + Mi + 2) 



the problem of the magnetic line blanketing can be simplified 
considerably by assuming that the magnetic field vector is ori- 
ented perpendicular to the line of sight, i.e. magnetic field is 
horizontal. In this case the effects produced by the polarized ra- 
diative transfer on the Stokes / are minimal and we can use the 
transfer equation for the non-polarized radiation, treating indi- 
vidual Zeeman components as independent lines. This approx- 
imation was u sed for the spe ctru m synthesis analysis of indi - 
vidual lines bv lTakedaldl993l) and lNielsen & Wahlgrenll2 002). 
The method is fully justified for weak lines (Stenflo 1994) 
and produces reasonable results for moderate and strong spec- 
tral features, although neglecting the polarized radiative trans- 
fer overestimates the magnetic intensification for the strongest 
spectral lines. At the same time, given the complex and di- 
verse angular dependence of the magnetic intensification of 
lines with different Zeeman patterns (iStift & Leo ne 2003), it is 
difficult to estimate the bias introduced by restricting our mod- 
elling to a horizontal field instead of dealing with a full range 
of the surface field orientations typical of magnetic CP stars. 

We modified the original line list by inserting additional 
spectral lines which correspond to individual Zeeman compo- 
nents of the anomalous splitting patterns. The wavelength shifts 
were calculated using Eq. (|5} and the oscillator strengths were 
established from the gf values of the original lines using the 
following formulas 

(gf)i„ = ~ S i P (gf)> 

\ (6) 

where the sum of the relative strengths is normalized to unity 
for each group of the Zeeman components 

N p N b N, 

£S,=J>,=J>, = 1. (7) 

!=1 !=1 1=1 

The last column in Table ^ shows the total number of the 
Zeeman components included in our line blanketing calcula- 
tions. 

3. Numerical results 

Using the LLModels code and treating the magnetic opacity 
as detailed above, we have calculated a set of model atmo- 
spheres with the effective temperature = 8000 K, 1 1 000 K, 
15 000 K, surface gravity logg = 4.0, metallicity [M/H] = 0.0, 
+0.5, +1.0 and magnetic field strength 0, 1,5, 10, 20 and 40 kG. 
For a few models with [M/H] - +1.0 we also considered 



the enhanced chromium abundance of [Cr] = +2.0 typical for 
some of the Cr-peculiar stars. The model atmosphere grid in- 
vestigated in this paper covers substantial part of the stellar pa- 
rameter space occupied by the magnetic Ap and Bp stars. 

The model atmosphere calculations were carried out on the 
Rosseland optical depth grid spanning from +2 to -6.875 in 
logT loss and subdivided into 72 layers. Convection was ne- 
glected, since it is generally believed that strong global mag- 
netic fields inhibit turbulent motions in the atmospheres of pe- 
culiar stars. Due the assumed absence of convection and be- 
cause of a direct inclusion of the magnetic intensification in 
the modelling of spectral lines, we adopted zero microturbu- 
lent velocity in all calculations with magnetic field. A possibil- 
ity to use pseudo-microturbulent velocity to mimic the effects 
of magnetic field on the line opacity was investigated by com- 
puting several model atmospheres with zero field and microtur- 
bulence in the range between 1 and 8 kms -1 . 

Numerical tests demonstrated that a sufficient accu- 
racy of the flux integration is achieved by using the 
following wavelength ranges for the spectrum synthesis: 
from 500 A to 50 000 A for r eff = 8000 K, from 500 A to 
30 000 A for r eff = 11 000 K and from 100 A to 30 000 A for 
Teff = 15 000 K. For all models we used a wavelength step of 
0. 1 A, which resulted in the total number of frequency points 
in the range between 295 000 and 495 000. 

At the first stage of creating models for each of t he stud- 
ied 7W , [M/ H] pairs we have used a standard atlas9 dKuruczl 
1993) model atmosphere for the spectral line preselection and 
as an initial guess of the model structure. The line selection 
threshold was set to 1 %. Experiments with the smaller selection 
cutoffs in the range from 10~ 4 % to 0.5% showed that, whereas 
including a larger number of lines results in a considerable in- 
crease of computing time, it leads to no more than 20-30 K 
difference of the temperature in the surface layers. 

Two different initial approximations were investigated for 
calculation of magnetic model atmospheres. The first one uses 
a model atmosphere with zero magnetic field, whereas in the 
second one a new model is calculated using a converged model 
with a weaker field strength (i.e. the OkG model is used as an 
initial guess for the 1 kG model, then the 5 kG model is cal- 
culated starting from the 1 kG model, etc.). We found that the 
results of the two approaches are identical, but the second one 
requires far less computing time due to a faster convergence. 
Consequently, this iterative approach was used for the calcula- 
tion of the magnetic model atmosphere grid presented here. 

In the following sections we study the influence of the mag- 
netic line blanketing on the model temperature and pressure 
structure and then investigate the effects on the common pho- 
tometric and spectroscopic observables: spectral energy distri- 
bution, photometric colors in the Stromgren, Geneva and Aa 
systems and profiles of the hydrogen Balmer lines. 

3.1. Model structure 

In Fig.^ we show the difference between the temperature of 
magnetic and non-magnetic models as a function of the optical 
depth. The magnetic line blanketing typically increases tem- 



O. Kochukhov et al.: Stellar model atmospheres with magnetic line blanketing 



5 



500 
400 
300 
200 
100 


-100 
500 
400 
300 
200 
100 


-100 
500 
400 
300 
200 
100 


-100 




T eff = 11000K, [M/H] = 0.0 



T eff = 15000K, [M/H] = 0.0 



^ \ '. 



-4 -3 -2 

log x 





Fig. 1. Difference of temperature between magnetic and non-magnetic model atmospheres for the effective temperatures 
r eff = 8000 K, 1 1 000 K, 15 000K and metallicities [M/H] = 0.0, +0.5, +1.0. 



perature in the line forming region. This distortion becomes 
more pronounced in a stronger field and higher metal abun- 
dance. The most persistent effect of including magnetic field in 
the line opacity calculation is the heating of atmosphere in the 
optical depth range from log T ross = -2 to 1 . This phenomenon 
appears due to the enhanced line blanketing resulting from the 
additional absorption in Zeeman components. The additional 
opacity leads to a redistribution of the absorbed energy back 
to the lower atmospheric layers and, thus, to heating of the 
log r ross -2 + 1 region. This effect is present for all the models 
shown in Fig.[Q 

In contrast, the behaviour of the upper layers depends 
strongly on the model parameters. For coolest stars considered 
here (T e g = 8000 K) we obtain a decrease of temperature in 
the upper layers. For the r e fr = 11 000 K models this cooling 
is limited to a small range of depths centred at log t 10SS « -3.5, 
whereas the magnetic models computed for r e (j = 15 000 K 
display a significant increase of temperature throughout the at- 
mosphere. This behaviour is explained by the shift from vi- 
sual to UV of the wavelength interval which plays the most 
important role for the atmospheric energy balance. For the hot- 
ter models the line density in UV becomes very high due to 
Zeeman splitting and, therefore, even relatively high layers be- 
come opaque and are able to absorb additional incoming radi- 
ation and contributes to blanketing effect. 



The pressure difference between magnetic and non- 
magnetic model atmospheres is shown in Fig.|2] The pressure 
stratification generally follows the tendencies of the tempera- 
ture distribution: deviations from the non-magnetic models be- 
come larger with higher metal abundance and stronger mag- 
netic field. The pressure tends to decrease when magnetic ef- 
fects on the line opacity are taken into account. This is a con- 
sequence of a non-linear response of the model structure to the 
modification of the line opacity and the respective temperature 
increase due to the anomalous blanketing. 

3.2. Energy distribution 

The spectral energy distributions for the three values of effec- 
tive temperature and metal overabundance of [M/H] = +0.5 
and [M/H] = +1.0 are presented in Fig. [3] The theoretical en- 
ergy distributions for magnetic stars are compared with the cal- 
culation for the reference non-magnetic models. 

Three main anomalies are evident in the computed flux dis- 
tribution of magnetic stars. The most conspicuous signature of 
the magnetically modified line blanketing is a flux deficiency in 
the ultraviolet spectral region and the respective flux excess in 
the visual. The magnitude of the UV deficiency is a moderate 
function of the effective temperature and clearly increases with 
the magnetic field strength and metallicity. 
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Fig. 2. Difference of pressure between magnetic and non-magnetic model atmospheres for the effective temperatures 
r eff = 8000 K, 1 1 000 K, 15 000K and metallicities [M/H] = 0.0, +0.5, +1.0. 



The presence of a magnetic field changes the stellar flux 
distribution in opposite direction in the visual and UV regions. 
At short wavelengths magnetic star appears to be cooler in 
comparison with a non-magnetic object with the same funda- 
mental parameters, whereas in the visual magnetic star mimics 
a hotter normal star. The "null wavelength" where flux remains 
constant progressively shifts to bluer wavelengths as the stellar 
effective temperature increases. 

Finally, we find that some of the theoretical flux distribu- 
tions of magnetic stars exhibit depression in the 5200 A region. 
This spectral feature has a well-known counterpart frequently 
observed in the spectra of peculiar stars ( Kupka et al. 2003). In 
our theoretical calculations the 5200 A depression is prominent 
at lower T e ff but becomes rather small for hotter models. The 
magnitude of the depression increases with the magnetic field 
intensity and metal content of the stellar atmosphere. 

3.3. Colors 

We studied the influence of the magnetic line blanketing on the 
photometric colors in the Stromgren uvbyH/3, Geneva and Aa 
systems. The dependence of various photometric parameters on 
magnetic field strength is shown in Figs.0]and|5] 

All colors we re calculated using modified computer codes 

I 1 J — — I 

by Kurucz ( 1993), which take into account transmission curves 
of individual photometric filters, mirror reflectivity and a pho- 



tomultiplier response function. In contrast to the Kurucz's pro- 
cedures (Relv ea & Kuruczlll978l) which are based on the low- 
resolution theoretical fluxes, our synthetic colors are computed 
from the energy distributions sampled every 0.1 A, so integra- 
tion errors are expected to be small. 

Th e peculiar parameter Aa was introduced by iMaitzenl 
to measure the strength of the 5200 A depression. Our 
calculations of the synthetic Aa photometry were carried out 
using the tabulated filter transmission curves and the photo- 
multiplie r response fu nctions as discussed in the recent paper 
by Kunk a et alJ (120031) . In addition, the peculiarity parameters 
Z and A(Vi -G) of the Geneva photometric systems were com- 
puted according to lCramer & Maederldl979|) . 

As one expects, the behaviour of the photometric indices is 
closely related to the flux redistribution between the visual and 
UV regions and the presence of several flux depressions, such 
the one at 5200 A. The shift of the "null wavelength" leads 
to alternation of the Balmer discontinuity amplitude as mea- 
sured by the c\ index of the uvbyH/3 system (see Fig.0}. We 
found that for all metallicities a relation between c\ and the 
magnetic field intensity is almost linear for T e g = 8000 K and 
T e s = 15 000 K, with positive and negative correlation respec- 
tively. On the other hand, for T e ff = 11 000 K the c\ index does 
not show any trend with the magnetic field strength because for 
this temperature the flux increases by a similar amount on both 
sides of the Balmer jump. 
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Fig. 3. Synthetic energy distributions from UV to near IR region for effective temperatures T e g = 8000 K, 1 1 000 K, 15 000 K and 
metallicities [M/H] = +0.5, [M/H] = +1.0. To enable visual comparison of the different curves the original LLModels stellar 
fluxes were convolved with a Gaussian profile with FWHM = 15 A. The insets show energy distributions in the 5200 A region. 



A relation between changes of all photometric indices and 
the intensity of magnetic field depends strongly on the stel- 
lar effective temperature. For low T e ff photometric changes are 
very pronounced, whereas for hotter magnetic stars modifica- 
tion of the uvbyH/3 photometric observables is fairly small (ex- 
cept ci). The Geneva photometric colors show a similar de- 
crease of the sensitivity to the magnetic effects with increasing 
T e g . This picture agrees with the general temperature behaviour 
of the theoretical flux distributions evident in Fig. [3] 

The Aa photometric system, characterizing the strength 
of the 5200 A depression, has long been considered as one 
of the most useful peculiarity indicators for stars with non- 
solar chemical c ompo s ition and strong magnetic field (e.g., 
see Kupka et al. 120031 120041 and references therein). Results 
of our calculations of the synthetic Aa parameter are summa- 



rized in Fig. |5] and confirm its sensitivity to the metal abun- 
dance and magnetic field strength. However, similar to the ef- 
fects observed for other medium and narrow-band photometric 
indicators, Aa is most influenced by magnetic field at lower 
r e ff, whereas the trend of Aa vs. the field strength quickly sat- 
urates for higher temperatures. For example, for = 8000 K 
and [M/H] = +1.0 the Aa index shows 0.043 mag differ- 
ence between the OkG and 20 kG field strengths, for T s s = 
1 1 000 K this difference decreases to 0.026 mag and it drops to 
0.010 mag for the models computed with T e g = 15 000 K. 

The prominent temperature sensitivity of the behaviour of 
the synthetic Aa indices arises because the absorption in the 
5200 A region is dominated by the Fei and low excitation 
Fe ii lines. These features become weaker with increasing T e g, 
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Fig. 4. Difference between photometric color indices of magnetic stars in comparison with the predictions for the non-magnetic 
reference models. Effect of magnetic field on the photometric parameters of the Stromgren system is shown for b—y (solid line), 
nt\ (dotted line), c\ (dashed line) and HB (dash-dotted line). The peculiarity indicators in the Geneva system, A(V] - G) and AZ, 
are shown with open circles and crosses respectively. 



hence the depression is stronger and more sensitive to the field 
strength in cooler magnetic stars. 

All photometric indices become more sensitive to the mag- 
netic effects with increasing metallicity. This is most clearly 
visible for the Aa index computed for T e $ = 8000 and 1 1 000 K 
models. The respective panels in Fig. [5] demonstrate that, for 
instance, the difference between Aa values for the models with 
r eff = 8000 K, [M/H] = 0.0 and [M/H] = +1.0 increases from 
the initial 0.010 mag for zero field to 0.03 1 mag for B = 40 kG. 
This effect is even more dramatic for T^g = 11 000 K where 
the difference Aa[M/H]=+\.o - Aa[M/H]=o.o grows from 0.016 to 
0.041 mag. 

An anomalous concentration of individual species with 
abundant spectral lines may have important implications for 
the Aa value. For the models with scaled solar abundance 
the 5200 A depression is produced mainly by the Fe lines. 
Contribution of the lines of other elements becomes significant 
only when the respective species are strongly overabundant. As 
illustrated in Fig. EI increasing the Cr overabundance by a fac- 
tor of 100 relative to the sun results in 0.016-0.023 mag growth 
of the Aa index for T eB = 8000 K, [M/H] = +1.0 model and the 
fiel d strength 5-2 kG. 

lHauckl ill 974 has demonstrated that the A(Vi - G) index 
of the Geneva system could be used as a peculiarity parame- 



ter for Ap stars. We found that the behaviour of A(Vi - G) is 
comparable to that of the Aa index. The maximum change of 
the A(Vi - G) parameter relative to the non-magnetic case is 
0.060 mag for B = 40 kG, T eff = 8000 K and [M/H] = +1.0. 

Th e X, Y, Z parameters introduced by Cramer & Maeder 
(1980) are linear combinations of the Geneva photometric in- 
dices. The Z value was suggested to be an indicator of ch emical 
peculiarity and the strength of the surface magnetic field (North 
1980). We found that the maximum change of Z relative to the 
non-magnetic case equals to -0.044 mag for T e g = 8000 K, 
[M/H] - +1.0 and B = 40 kG. However, at higher tempera- 
tures Z becomes less sensitive to the magnetic field effects com- 
pared to A(Vi - G) and Aa. 

Finally, we note that none of the photometric indicators of 
magnetic stars proposed in the literature shows a linear trend 
over the whole range of the considered magnetic field strength. 
Relations between changes of the photometric indices and the 
magnetic field intensity show saturation effects for B > lOkG 
This saturation can be offset to stronger fields by increasing the 
number of spectral lines affected by the Zeeman splitting and 
intensification, which explains greater magnetic sensitivity of 
the synthetic photometry computed for the models with lower 
r e ff and higher metal abundance. 
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Fig. 5. The magnitude of the Aa photometric index as a function of the magnetic field strength for T e g = 8000, 1 1 000 and 
15 000 K. Different curves show calculations with metal abundances [M/H] = 0.0 (open circles), [M/H] = +0.5 (triangles) and 
[M/H] = + 1.0 (filled circles). Open squares illustrate the effect of increasing Cr overabundance to [Cr] = +2.0inther e ff = 8000, 
[M/H] = +1.0 models. 



3.4. Hydrogen line profiles and metal lines 

Profiles of the Ha, HB and Hy hydrogen Balmer lines were cal- 
culated using the Synth program (Piskunov 1992). This spec- 
trum synthesis code incorpo rates the new hydrogen line Stark 
broadening computations bv lStehlel JT994). as well as the im- 
proved treatm ent of the hydr ogen line self -broadening devel- 
oped by Barklem et al 1 (120001) . The comparison between the 
magnetic and non-magnetic profiles of the H/3 line calculated 
for models with 10 times the solar metal abundance is pre- 
sented in Fig.|6] 

We found that the changes in the atmospheric structure due 
to the magnetic line blanketing do not have a strong influence 
on the hydrogen line profiles. For HB the maximum change 
relative to the non-magnetic model amounts to about 3% of the 
continuum level for B = 40 kG, but it does not exceed 1% for 
Z? < 10 kG, which is a more typical surface field strength for 
the majority of magnetic CP stars. 

To estimate an effect of the anomalous atmospheric struc- 
ture of the magnetic models on the profiles of individual metal 
lines we computed synthetic spectra in the 4000-6000 A re- 
gion for the zero field and the 10 kG field models and [M/H] = 
+ 1.0. Since at this point we are interested to study only the 
effects of the model structure, the spectrum synthesis calcula- 
tions were carried out disregarding magnetic field. We found 
that the depth of absorption lines is modified in the spectra 
computed with the magnetic models. The strongest effect is 
found for T e ff = 8000 K, where in the presence of the 10 kG 
field many medium strength and weak lines become sytemati- 
cally shallower by 4-6% of the continuum level relative to the 
calculation with the zero field model. This discrepancy is re- 
duced to less than 2-3% for the T eff = 11 000 and 15 000K 
models. For these hotter atmospheres the difference between 
the line profiles computed for the magnetic and non-magnetic 
models can be both negative and positive. For many lines the 
main discrepancy is concentrated in the line wings, with the 
magnetic models predicting slightly narrower line profiles. 



3.5. Stellar atmospheric parameters and bolometric 
correction 

The overall modification of the energy distribution of magnetic 
stars as well as the presence of localized absorption features 
may bias photometric determination of the stellar atmospheric 
parameters. As a part of the investigation in this paper, we ver- 
ified an influence of the magnetic line blanketing on the photo- 
metric determination of T e ff and log g. 

The TempLogg code dRogersll995l) was applied to the syn- 
thetic uvbyHB colors calculated as described in Sect. 13.31 We 
found that for = 8000 K the maximum difference between 
r e ff and logg obtained for non-magnetic models and the corre- 
sponding parameters for the B = 40 kG models is -190K and 
-0.27 dex respectively. For the r e ff = 11 000 K models the tem- 
perature discrepancy is maximal (+86 K) for normal composi- 
tion and decreases with increasing metal abundance, whereas 
the difference of log g peaks at -0. 1 1 dex. For the hottest mod- 
els (r eff = 15 000 K) the r eff difference is +541 K, but logg is 
not affected. It is clear that, even in the case of extreme mag- 
netic field strengths, the bias in determination of the fundamen- 
tal parameters using the Stromgren photometry does not result 
in a systematic deviation of T e g and log g beyond the usual er- 
ror bars assigned to the photometrically determined stellar pa- 
rameters. 

Somewhat different results emerged from our analysis of 
the stellar parameters determined using the Geneva photomet- 
ric sys tem. Effective temp erature was derived applying the 
code of lKiinzli et alJ (119971) to the synthetic Geneva colors and 
was found to increase systematically with the field strength for 
all three considered T e s values. The bias introduced by mag- 
netic field reaches 100-200 K for T eS = 8000-11 000 K and 
900 K for T e ff = 15 000 K. The surface gravity is modified by 
up to 0.2 dex. Thus, a higher sensitivity to the metallic line 
absorption makes the standard Geneva photometric calibration 
a less robust method for determination of the atmospheric pa- 
rameters of magnetic stars compared to the calibration in terms 
of the uvbyHB colors. 
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Fig. 6. Comparison between the synthetic HB profiles computed for [M/H] — +1.0 and different values of the magnetic field 
strength and r e g\ The lower panels show normalized profiles, whereas the upper plots illustrate the difference between HB 
calculated for the models with substantial magnetic line blanketing and the reference non-magnetic model. 



An increase of the line blanketing due to magnetic effects 
and overabundance of metals lead to a redistribution of the stel- 
lar flux from UV to visual. Consequently, magnetic CP stars 
appear brighter in V compared to normal stars with the same 
fundamental parameters. Using the theoretical flux distribu- 
tions calculated with the LLModels code we studied the ef- 
fect of magnetic field and metallicity on the bolometric correc- 
tion (BC). For all stellar models this parameter decreases rel- 
ative to the non-magnetic case by up to 0.054-0.089 mag for 
B = 40 kG and by no more than 0.025 mag for B < 10 kG. The 
BC for the r e ff = 15 000 K models is the most sensitive to mag- 
netic field. At the same time, even in the absence of the field, an 
increase of metal abundance leads to a substantial modification 
of BC. It decreases by 0.079-0. 115 mag when metal abundance 
is changed from solar to [M/H] = + 1 .0. 

3.6. Enhanced microturbulence models 

Until the present investigation the most common approach 
to account for the magnetic line blanketing in the model at- 
mosphere analyses of magnetic CP stars w as to use an in- 
creased value of t he micro turbulent velocity (Muth samll 19791 
Kup ka etafll2004l) . For the magnetic field B such a pseudomi- 
croturbulent velocity is computed according to the following 
equation (e.g., Kupka et al. 1996): 

Snag =4.66xl0- I3 C^eff.B, (8) 

where c is the speed of light in kms" 1 , the field strength B 
is measured in gauss, A is the wavelength in A and g e g is the 



effective Lande factor. In previous model atmosphere studies 
A * 5000 A and g e ff = 1.0—1.2 were typically adopted, which 
gives £ m ag ~ 4 and 8 km s" 1 for the magnetic field strength of 5 
and 10 kG respectively. 

Having introduced a more realistic implementation of the 
magnetic intensification in modelling of stellar atmospheres, 
we are in the position to verify the extent to which models 
with the magnetic pseudomicroturbulence are able to match the 
properties of our more sophisticated but, admittedly, also more 
computationally expensive models. To study the performance 
of the enhanced microturbulence models, the model structure, 
energy distribution and synthetic Aa indices computed for zero 
field and £ mag = 1-8 km s" 1 were compared with the respective 
parameters for the B = 5 and 10 kG models from our main grid. 

It appears that the enhanced microturbulence is indeed able 
to partially mimic the behaviour of magnetic models, such as 
the energy redistribution from UV to visual and the heating 
in certain atmospheric layers. However, an exact quantitative 
match requires adopting a different £ mag parameter for each 
considered quantity. For instance, to reproduce the tempera- 
ture structure and energy distribution of the T e g = 8000 K 
models one has to use (depending on metallicity) £ mag = 2.1- 
2.6 kms" 1 and 3.7-4.1 kms" 1 for B = 5 and 10 kG respec- 
tively, which is roughly a factor of two smaller microturbu- 
lence than predicted by Eq. At the same time, matching 
the anomalous absorption in the 5200 A region requires a 2- 
3 kms" 1 higher value of £ mag . Furthermore, models' sensitiv- 
ity to the magnetic line blanketing decreases with T e g faster 
than the influence of microturbulence. Consequently, for the 
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T e g = 15 000 K models one has to adopt nearly twice as weaker 
£mag to simulate effect of the same field as in cooler models. 

Results of the calculations in this section demonstrate that a 
properly chosen value of £ mag can only be used as a very rough 
guess of the effects due to magnetic field. However, it is evident 
that the models with enhanced microturbulence are not suitable 
for detailed analysis of magnetic stars and cannot be used for 
modelling individual features of the stellar flux distributions. 

4. Conclusions and discussion 

Construction of the model atmospheres of CP stars is consid- 
erably more complicated in comparison to that of normal stars 
because such modelling has to take into account and explain a 
range of phenomena not present in the atmospheres of normal 
stars. The main and most important distinction of magnetic CP 
stars from the other Main Sequence objects is the presence of a 
strong global magnetic field. It plays a central role in producing 
chemical abundance anomalies and inhomogeneities and leads 
to a peculiar atmospheric structure, unusual line profiles and an 
anomalous flux distribution. 

The Zeeman splitting of spectral lines is one of the principal 
effects taking place in the magnetic stellar atmosphere. In this 
study we have computed a grid of model atmospheres of pecu- 
liar A and B stars accounting for the magnetic line blanketing. 
Our model grid has covered the range of possible atmospheric 
parameters of CP stars, which allowed us to analyze the be- 
haviour of the model structure and observed characteristics for 
different magnetic field, abundances and effective temperature. 
We have used the new model atmosphere code LLModels and 
up-to-date compilations of the atomic line data from the VALD 
database to insure an accurate modelling. 

A somewhat simplified "horizontal" model of magnetic 
field allowed us to reduce computational costs by solving the 
radiative transfer equation for the non-polarized radiation. This 
approach does not fully reflect the influence of magnetic field 
on the formation of spectral lines but, nevertheless, it represents 
the most accurate treatment of the Zeeman effect in model at- 
mosphere calculation used so far. 

Our models show that the enhanced line blanketing due to 
the magnetic intensification of spectral lines produces heating 
of the certain atmospheric layers and this effect increases with 
effective temperature, magnetic field strength and metal abun- 
dance. 

We found that the model atmospheres with the magnetic 
line blanketing produce fluxes that are deficient in the UV re- 
gion by 0.2 - 0.3 mag for B = lOkG and about 0.3 - 0.8 mag 
for B = 40 kG. Moreover, the presence of a magnetic field leads 
to the flux redistribution from the UV to visual region due 
to backwarming. This property of the theoretical models is in 
agreement wi th those observations of CP stars llLeckr one 1974; 
iMolnarl 1 973H.TamaJ 1 977l) which have revealed the presence of 
the flux redistribution and demonstrated that the flux changes 
in the visible and ultraviolet regions are inversely correlated. 

The energy distributions of the calculated model atmo- 
spheres show that an increase of effective temperature shifts 
the "null wavelength" where flux remains unchanged to shorter 
wavelengths. The results of our modelling suggest that the "null 



wavelength" approximately falls at 3800 A for J e ff = 8000 K, 
at 2900 A for r eff = 1 1 000 K and at 2600 A for T eB = 1 5 000 K 
and slightly depends on metallicity and field strength. This fact 
is in concordance with t he obs ervational data and numerical ex- 
periments by iLeckrone et alJ II 19741) who examined enhanced 
opacity models and found that the "null wavelength" shifts to- 
ward UV when effective temperature is increased. 

Magnetic CP stars demonstrate several flux depressions in 
the visual and ultraviolet region unlike other stars. The most 
prominent feature in the visual is the flux depression centred at 
5200 A. Our numerical results show that the strong feature near 
5200 A appears for lower effective temperature and vanishes 
for higher effective temperature and its amplitude increases 
with the magnetic field strength and metallicity. Comparison of 
the synthetic Aa indices with typical values of this parameter 
observed for cool magnetic stars demonstrate that our calcula- 
tion are able to reproduce the whole range of the 5200 A de- 
pression strength. Furthermore, our results argue that for mod- 
elling the 5200 A feature in cool CP stars an accurate treatment 

of magnetic field is of no less im portance compared to using 

j II I 

individual abundances (Kupka et al. 2004). A 10 kG magnetic 

field changes the value of the Aa parameter by as much as 20- 

30 mmag, which is substantial even compared to the largest 

observed Aa values. 

Since, historically, photometric data are widely used for 
the stellar classification and atmospheric parameter determi- 
nation, we have calculated indices in the uvbyHB and Geneva 
photometric systems as well as the peculiar indices Aa, Z and 
A( Vi - G). We found that changes of most photometric param- 
eters due to the influence of magnetic field are noticeable for 
low effective temperatures, whereas for hotter stars sensitivity 
to magnetic field is reduced considerably. This behaviour of the 
theoretical models a grees with the flux distribution properties 
and with the study bv lCramer & Maeden I I 19801) . who noted that 
the effects of magnetic field on the energy distribution become 
less important with an increase of temperature. 

The peculiar indices Aa, Z and A(Vi - G) are frequently 
used for identification and assessment of the properties (in par- 
ticular the surface field strength) of CP stars. We examined the 
influence of magnetic field on these indices. We found that the 
magnetic modification of the peculiar photometric parameters 
is relatively small (maximum change of Aa is 0.065 mag for 
40 kG field) and the relation between their values and the mag- 
netic field strength is not linear due to the saturation effects for 
stronger field (B > 10 kG). In any case, the Aa parameter ap- 
pears to be the most sensitiv e to magnetic field, where as Z and 
A( Vi - G) are less sensitive. Cramer & Maeden Jl980h pointed 
that the saturation of the photometric effects for large fields is 
an observational fact and has to be explained by a model atmo- 
sphere analysis. Our numerical results reproduce and provide a 
clear explanation of the behaviour of the photometric indices. 

Finally, we investigated the question to what extent the 
anomalous magnetic opacity affects the photometric determi- 
nation of the atmosphere parameters of CP stars. We found that, 
although an enhanced metal abundance and magnetically mod- 
ified line opacity produce some changes of the model atmo- 
sphere structure and the flux distribution, these changes are not 
very large and do not result in a strong modification of the op- 
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tical photometric colors. Consequently, the model atmosphere 
parameters derived using the photometric calibrations for nor- 
mal stars are not far from their true values. Thus, despite the 
presence of several local anomalies of the stellar energy distri- 
bution and the overall flux redistribution from the UV to visual 
region, peculiar stellar models do not appear to mimic standard 
models with a systematically different T e g. Standard calibra- 
tions of the Stromgren photometry appear to be especially ro- 
bust for the determination of the parameters of magnetic stars. 
Thus, we conclude that by itself the magnetic opacity does not 
introduce significant errors in the photometric estimates of CP- 
star parameters. 

We also showed that the hydrogen Balmer line and absorp- 
tion features of metals are not very sensitive to the changes in- 
duced by the magnetic line blanketing on the model atmosphere 
structure. Although the effect on the overall flux distribution 
may be noticeable, the changes in the depths of metal lines do 
not exceed a few % of the continuum level for the 10 kG field, 
which justifies neglecting magnetic line blanketing in the rou- 
tine abundance and spectrum synthesis studies of the majority 
of magnetic CP stars. 
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